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ABSTRACT 

Quantum  lattice  gas  algorithms  are  developed  for  the  coupled-nonlinear  Schrodinger  (coupled-NLS)  equations, 
equations  that  describe  the  propagation  of  pulses  in  birefringent  fibers.  When  the  cross-phase  modulation  factor  is 
unity,  the  coupled-NLS  reduce  to  the  Manakov  equations.  The  quantum  lattice  gas  algorithm  yields  vector  solitons  for 
the  fully  integrable  Manakov  system  that  are  in  excellent  agreement  with  exact  results.  Simulations  are  also  presented 
for  the  interaction  between  a  turbulent  2-soliton  mode  and  a  simple  NLS  2-soliton  mode.  The  quantum  algorithm 
requires  4  qubits  for  each  spatial  node,  with  quantum  entanglement  required  only  between  pairs  of  qubits  through  a 
unitary  collision  operator.  The  coupling  between  the  qubits  is  achieved  through  a  local  phase  change  in  the  absolute 
value  of  the  paired  qubit  wave  functions.  On  symmetrizing  the  unitary  streaming  operators,  the  resulting  quantum 
algorithm,  which  is  unconditionally  stable,  is  accurate  to  0(e2) . 

Keywords  :  vector  solitons,  coupled-NLS,  Manakov  equations,  quantum  lattice  gas  algorithms,  turbulence 

1.  INTRODUCTION 


One  of  the  most  fundamental  equations  of  soliton  physics  is  the  nonlinear  Schrodinger  equation  (NLS) 


+  e[M]  V  =  0 


0) 


which  can  describe  the  propagation  of  optical  solitons  in  an  ideal  fiber.  However,  in  reality,  optical  fibers  are 
birefringent  and  a  single  mode  fiber  can  support  two  distinct  orthogonal  polarizations.  One  of  these  modes  is  the  O- 
mode  and  has  constant  refractive  index  along  the  ray  path  while  the  other  mode  is  the  X-mode  whose  refractive  index 
varies  along  the  ray  path.  Thus  the  pulses  will  travel  down  with  slightly  different  speeds  along  the  two  orthogonal 
polarization  axes.  Proceeding  from  the  electromagnetic  wave  equation  for  the  electric  field  E  in  a  dielectric  medium 
1  92E  32P 

V  E  -  TTY  =  “Mo  TT  (2) 

c  dt  dt 

where  the  polarization  P  represents  the  linear  and  nonlinear  induced  polarization,  it  can  be  shown  [1]  that  the  slowly 
varying  amplitudes  of  these  modes  satisfy  the  coupled-NLS  equations 
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Here  g  is  a  positive  parameter  and  B  is  the  cross-phase  modulational  coefficient,  B =  (2  +  2 sin2  6>) /(2  +cos2#), 


where  0<9<n/2  is  the  birefringence  ellipticity  angle  [2].  The  coupled-NLS  equations  (3)  are  generally  nonintegrable, 
except  for  the  special  case  of  £=1  in  which  they  reduced  to  the  completely  integrable  Manakov  equations  [3].  The 
Manakov  equations  are  completely  integrable,  with  soliton  solutions.  For  a  linearly  birefringent  fiber  5=2/3.  The 
coupled-NLS  equations  also  arise  in  spatial  beam  propagation  in  crystals  [4],  and  under  special  circumstances  5  =  1, 
yielding  the  soliton  Manakov  system..  These  Manakov  equations  also  describe  propagation  in  randomly  birefringent 
fibers  [5]. 


Here  we  develop  a  quantum  lattice  gas  representation  of  the  more  general  counled-NT  s 
coupling  potentials  V1  and  V2  : 
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which  is  unconditionally  stable  and  ideal  for  implementation  (and  parallelization)  on  a  hybrid  quantum-classical 
computer  [6]  as  well  as  on  a  classical  computer.  This  representation  is  a  generalization  of  our  earlier  work  [7]  on 
quantum  algorithms  for  the  simple  NLS  and  KdV  solitons,  which  was  based  on  the  quantum  algorithm  for  the  (linear) 
Schrodinger  equation  [8],  These  quantum  lattice  gas  algorithms  can,  in  principle,  be  modeled  on  a  liquid  NMR  quantum 
computer  [9-13].  The  exponential  speed-up  over  classical  computers  arises  from  the  quantum  entanglement  of  qubits.  In 
Sec.  2  we  briefly  summarize  the  quantum  lattice  gas  algorithm  for  4  qubits/lattice  The  unitary  collision  and  streaming 
operators  entangle  just  pairs  of  qubits.  The  local  coupling  of  all  the  qubits  is  achieved  by  introduce  a  phase  change  into 
the  wave  functions.  This  phase  change  (for  each  qubit  pair)  allows  for  the  introduction  of  an  "external"  potential  and  the 
modeling  of  vector  soliton  physics.  For  the  liquid  NMR  computer,  one  needs  only  hold  the  phase  coherence  for  16  gate 
operations  before  measurements  are  needed  in  order  to  determine  the  potentials.  Vector  soliton  collisions  are  then 
studied  with  this  algorithm  in  Sec.  3,  for  both  the  Manakov  system  as  well  as  for  the  nonintegrable  case  of  1  and  for 
vector  soliton  turbulence.  Excellent  agreement  is  found  for  both  the  vector  1 -soliton  and  2-soliton  solutions  of  the 
Manakov  system. 

2.  QUANTUM  LATTICE  GAS  REPRESENTATION  FOR  COUPLED-NLS 

We  first  discretize  the  spatial  domain  into  L  lattice  nodes  so  that  the  wave  function  i jf  has  a  basis  ket  |  x1 }-  expansion 

L 

k)  =  I>i|*i)  ,  with  c1=(x1[i^)  (5) 

1=1 

being  the  probability  amplitude  associated  with  the  ket  |  xx  ).  For  the  coupled-NLS  equations  we  introduce  4  qubits  at 
each  lattice  node  1 :  ,  a  =  0,1, 2, 3.  Each  qubit  is  a  two-level  quantum  system 

|^a)  =  aa  |0)+  /?a|l)  ,  with  ,  a  =  0...3  ;  1=1. ...L  (6) 

In  the  number  representation  one  can  employ  binary  indexing  for  the  basis  set 

\nl0n\  n‘n‘  ....  n\nLxn\n\)  where  «(u=O0rl  >  V  3,  and  n23=0  or  1  ,  Vl  (7) 

which  factors  since  only  the  qubits  |<7o,i)  w'll  be  entangled  with  each  other  and  similarly  for  the  entanglement  of  the 
qubits  ^2,3}  at  each  spatial  site: 


1  1 

L  L\ 

1  1 

L  L\ 

nQnx  . 

-  »o  n\  ) 

n2n3  . 

••  *2  n3  ) 

(8) 

We  can  restrict  ourselves  to  the  one-particle  sectors  for  qubits  |<7o,i)  ant^  qubits  |  ^23)  *n  which  the  basis  set  elements  in 
Eq.  (8)  have  only  one  Bqj  that  is  1  while  all  the  other  rigg  are  zero,  and  similarly  for  n2,3-  In  essence,  for  the  two  qubits 
|^o, i)>  there  are  2L  such  elements  which  can  be  labeled  using  the  binary  scheme  1 221+a^,  1=1. ..L  ,  a=0,l  so  that  in 

this  binary  scheme  the  wave  function  corresponding  to  the  two  qubits  1 9o,i) 

H,  -  X  X  {*„  b”-1*')  ■  (8) 

1=1 0=0 

Interference  effects  between  the  two  qubits  |<?o,i)  at  ea°h  position  ket  |xx)  can  be  achieved  by  setting 

C1  =l32l-l  +  £21  (9) 

Similarly  for  the  interference  effects  between  the  two  qubits  |  ^23)-  To  recover  the  coupled-NLS  equations  we  do  not 
need  to  enforce  interference  effects  between  all  four  qubits  at  the  spatial  node  Ix-J.  Hence,  for  simplicity,  we  now  just 
concentrate  on  the  application  of  the  unitary  collision  and  streaming  operator  to  the  two  qubits  j  g0i)  —  exactly  the  same 
unitary  and  streaming  operators  will  be  applied  to  the  two  qubits  |  #2  3  ^ . 
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2.1  Unitary  collision  operator  for  the  qubit-pair 

To  evolve  the  wave  function  for  the  two  qubits  |cy0  J  in  time,  a  quantum  unitary  operator  C  is  constructed  from  a  tensor 
product  of  quantum  gates,  each  independently  applied  on  a  site-to-site  basis 

_  L  ^ 

C=®Ui  (10) 

1=1 

Quantum  entanglement  arises  from  the  on-site  local  unitary  collision  operator  U i  acting  on  the  2  qubits  |f/0l^  at  each 

node;  i.e.,  U  i  acts  on  the  4  on-site  basis  kets  for  j<y0li)  : 

|0)®|0)=|l000),  |0)®|l)s|0100)  ,|1)®|0)=|0010)  ,  |1)®|1)=|0001) 

In  particular,  U i  acts  on  the  on-site  ket  |v)  =  |0)®  |l)  +  |l)®  |0)  =|0110).  A  local  equilibrium  can  be  associated  with 
this  on-site  unitary  collision  operator  if  |v)  is  an  eigenvector  of  U  i  with  unit  eigenvalue  :  U  i  |v)  =|  v) . 

To  recover  NLS,  we  introduce  the  square-root-of-swap  gate  Ui=U  on  a  site-by-site  basis  [8] 

ho  o  ol 


n  1-i  1+j  „ 

0  —  —  0 
2  2 
1+j  1 - i 
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In  the  number  representation  UNLS  acts  on  the  kets  |221  and  |221^  for  the  two  qubits  |ry01^  on  the  site  {x^ . 

Moreover,  a  Hamiltonian  representation  for  this  unitary  quantum  gate  UNLS  can  be  achieved  from  appropriate  tensor 

( 0  n  (0  -A  fl  O'} 

products  ofthe  Pauli  spin  matrices  <TX  =1  q  I  ,  Oy  - 1  .  I  ,  and  az  =1  ^  forqubits  ’0’  and  ’1’ : 

Unls  =exP  Y  exP  ~J +  a'j*l  +  azaz)]  (]2) 

with  U^LS  =  / ,  the  identity  operator  so  that  U^LS  |v)  =  |v). 

2.2  Unitary  Streaming  Operator  for  the  qubit-pair  |  g'o.i ) 

The  next  step  of  the  quantum  lattice  gas  algorithm  is  to  stream  the  post-collision  on-site  ket  for  qubits  J  ^o,i)  t0  nearest 

neighbor  sites.  The  (unitary)  streaming  operator  Si  is  defined  as  a  global  shift  to  the  right  of  only  qubit  |<?0)on  each 
lattice  node,  i.e., 

Si =n%2m  03) 

1=1 

where  X 21-1 21+1  's  independent  of  1  and  in  the  number  representation  shifts  the  amplitude  of  the  ket  |221-1^  on  site 
{xx}  to  the  ket  J 2?1+1  j  on  site  {x1+1}for  the  qubit-pair  |<y0,i)  •  In  matrix  form,  the  qubit  streaming  operator  %  is  a 

22  X  22  (unitary)  permutation  matrix 
'1  0  0  0N 

0  0  10 

X  =  •  (14) 

A  0  1  0  0  v  2 

k0  0  0  lj 

Thus,  for  the  two-qubit  pair  |?o,i)>  U  operates  on  the  on-site  qubits  while  %  operates  on  the  1st  qubit  of  neighboring 
sites.  Hence  the  total  collision  matrix  C  does  not  commute  with  the  streaming  operator  Si.  To  symmetrize  the 
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algorithm,  we  also  introduce  the  streaming  operator  S->  that  performs  a  global  shift  to  the  right  of  the  2nd  qubit  on  all  the 
lattice  nodes  for  the  qubit-pair  |qu  ; 

L-l  ^ 

■^2  =  IT  £21.21+2  (15) 

1=1 

2.3  Introduction  of  the  Potential  Field  for  the  qubit-pair  |  q0 , ,  'j 

It  has  been  shown  [14,  8]  that  the  effect  of  an  external  potential  ^(x)  can  be  modeled  by  the  introduction  of  a  local 
phase  change  to  the  system  wave  function  for  the  qubit-pair  |  <?o,i)-  Let  Q\  denote  the  wave  function  for  the  qubit- 
pair  |  q q, ] ) ,  and  Q2  for  the  qubit  pair  |  <72,3)-  Thus 

Qi{x,t)  ->  exp[  i !/[(*)  At]  Qi(x,t)  (16) 

where  A  t  is  the  time  advancement  after  each  step  of  the  algorithm. 

Similarly,  for  the  qubit-pair  |  <?2,3 ^ -  One  can  introduce  another  external  potential  V2(x)  by  a  local  phase  change  in  the 
system  wave  function  Q2  for  the  qubit-pair  ^2,3)  : 

Q2(x,t)  -1  exp[ i  V2(x)  At]  Q2(x,t)  (17) 

2.4  Quantum  algorithm  for  coupled-NLS 

We  simultaneously  apply  the  collide-stream  sequence  of  unitary  operators  to  the  qubit  pairs  |  q0  Ij  and  |  q2  3j 

foAt  +  At)^  ffs2 c  s2c  sic  s2c  .sic  5,C  Sf c  S,c]  e,(t)N 

.  A  =  .  ~  _ _  (18) 

KQiV  +  At))  k'c  S2C  S2C  S2C  .  Si  C  5,C  S\  C  S,  C  Q2{t) 

VL  / 

where  Sf  is  the  transpose  of  St,  with  Sj  5,  =/,  i=l,2  and  C  is  based  on  the  unitary  collision  operator  UNLS,  Eq.  (11). 

We  then  apply  the  phase  change  transformations  (16),  (17)  following  the  collide-stream  sequence  : 

(QA  ^  Wv>(l  2,1,122 1 )  Atl Qi  j  (19) 

le2J  [exp^de,  |, | e2 1 ) At], e2J 

where  the  potential  fields  are  required  to  be  functions  of  the  wave  functions  themselves  V,  =  vTI  Qx  |,|  Q2  |) 


The  continuum  limit  is  defined  by  scaling  the  spatial  shift  between  neighboring  nodes  to  be  0(e),  the  time  advancement 
A r=(?(e2]  and  the  potentials  -e2  V)[|  Qx  |,| Q2  |] .  In  the  limit  £  — >  0,  it  can  be  shown  using  Mathematica  that  the  Eq. 
(18)-(19)  sequence  yields  the  coupled-NLS  equations 

l~dt~  +  ~d^  +  ’l^2l]  •  2i  =  0  +  0(e2)  , 

^  +  V2[|Q2Nef|]  •  fi2  =  0  +  o(e2)  ,  as  e—>  0  (20) 

with  error  of  O^e2).  Equation  (20)  holds  for  any  choice  of  potentials  Vj  and  V2. 

It  is  now  clear  how  to  extend  this  analysis  to  a  system  of  N-coupled  NLS  equations.  N-coupled  NLS  equations  arise  in 
the  study  of  beam  propagation  in  a  Kerr-like  photorefractive  medium,  which  typically  exhibits  very  strong  nonlinear 
effects  with  extremely  low  optical  powers.  Partially  coherent  solitons  have  been  observed  through  excitation  by  partially 
coherent  light  [15  ]  as  well  as  by  ordinary  incandescent  light  bulb  [16). 
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3.  QUANTUM  LATTICE  GAS  SIMULATIONS  FOR  COUPLED-NLS 
3.1  Vector  soliton  Manakov  system 

We  first  test  our  quantum  algorithm  on  the  fully  integrable  Manakov  system 

^  +  2/t[| 2,|2  +  \q2\2].  q,  =  o  ,  ^  +  2/t[|e2|2  +  |e,|2].  e2  -  o  (21) 

with  the  cross-phase  modulation  coefficient  B  =  1.  The  I -soliton  vector  solution  is  [17] 

cc  R  ~  R  B 

2iU.f)  =  — exp -—  +  it],  SechriR+—  ,  Q2(x,t)  =  ^Qx(x,t)  (22) 

where  T]=k(x  +  ikt)  and  a, (3 ,k  are  arbitrary  complex  parameters  with  subscripts  R  and  /  denoting  the  real  and 
imaginary  parts.  R  =  ln[u(ja  |2  +1  (3 12)  /  (k  +  k  *)2|  is  real  and  kR  ^  0. 


Q  (x  ,  t) 


1000  2000  3000  4000  5000 
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Fig.  1  The  Vector  2-soliton  Manakov  solution  [B  =  I]  for  (a)  pre-collision  and  (b)  post-collision  times  for  the  two  orthogonal 
polarization  direction  0|(.v.;)  and  Q^(x.l).  Initially  (  t  =  r(, ).  the  2-solitons  of  the  Qt -mode  are  at  x  =  61 5  and  x  =  4605  with  the 

larger  amplitude  soliton  moving  — »  while  the  lower  amplitude  soliton  has  the  same  speed  but  moving  For  the  2-solitons  of  the 
Q2  -  mode,  the  smaller  amplitude  soliton  at  x  =  615  moves  in  — >  while  the  larger  amplitude  sgliton  at  x  =  4605  moves  in  <— .  The 
2-vector  soliton  collision  is  elastic  as  seen  at  the  snapshots:  t0  =  0,  t{  =5.5  K  ,  t2  =  10  K  (collision  -  strong  wave  function 
overlaps),  r3  =  12.5  K  and  r4  =  20  K  iterations. 


For  non-overlapping  solitons,  the  initial  ( tQ)  2-soliton  states  for  the  Qx-  and  (22-modes  reduce  to  the  simple  1-soliton 
wave  functions  of  the  form  given  by  Eq.  (22).  The  vector  elastic  collision  of  these  2-soliton  modes  is  clear. 


3.2  Integrable-Nonintegrable  Coupled-NLS 

We  now  consider  the  following  evolution  equations  for  the  Ql  -  and  Q2  -modes: 
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with  the  same  initial  conditions  as  in  the  Manakov  system: 
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In  our  earlier  work  [7],  for  the  scalar  nonintegrable  NLS  for  the  uncoupled  Q2  -mode,  the  Q2  -mode  exhibited  soliton 
turbulence  [18].  Here  we  shall  see  that  the  coupling  of  these  two  modes  will  drive  the  Ql  -mode  out  of  its  quasi-soliton 
state  because  of  its  coupling  to  the  turbulent  Q2  -mode. 
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Fig.  2  The  initial  pulses  for  the  Q,  -  and  Q2  -  modes  ( t0  =  0)  which  are  the  2-soliton  vector  solutions  for  the  Manakov  system, 
Fig.  1.  At  =  0.7  Jf,  the  Qt  -mode  retains  the  vector  2-soliton  structure  while  the  Q2  -  wave  function  undergoes  very  significant 
contraction,  although  the  group  velocities  of  these  pulses  are  unaffected  by  the  nonintegrable  coupling  potentials  and  remain  equal  to 
the  Manakov  group  velocities. 

The  quasi-integrable  £?i  -mode,  in  the  early  time  evolution,  retains  its  Manakov-like  soliton  solution  (Fig.  2)  for  early 
times  ( t0+  =  0.7  K)  because  of  the  weak  coupling  to  the  nonintegrable  Q2  -  mode  which  rapidly  loose  their  Manakov- 
like  structure  and  become  more  localized  structures.  However,  the  structures  retain  their  Manakov-like  group  velocities. 
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While  the  nonintegrable  Q-,  -  mode  retains  its  sharp  structures  (Fig.  3  with  q  =  1.4  K).  but  small  amplitude  sidebands 
form  and  then  start  to  rapidly  move  away  from  the  sharp  structure  ( t,  -5.5 K).  The  coupling  of  the  (2,-  and 
Q2  -modes  now  distorts  the  quasi-Manakov  (2|  -solitons  with  slight  contraction  and  higher  peak  together  with  the 
emission  of  very  weak  sidebands  (Fig.  3  ). 


Fig.  3  The  evolution  of  |g,  |  ,  \Q2 1  at  times  tl  =  \AK  and  t2  =5.5 K .  The  quasi-Manakov  Qy  -  solitons  are  slightly  compressed 
due  to  the  coupling  with  the  nonintegrable  Q2  -structures  which  become  very  sharp  and  emit  weak  sidebands. 

As  these  mode  develop  in  time,  Fig.  4,  one  notes  that  the  all  the  structures  in  the  Q2  -mode,  both  the  very  sharp  large 
amplitude  peaks  and  the  broadband  ’noise’  with  rapidly  moving  weaker  fine  structures,  behave  like  quasi-solitons  :  i.e., 
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Fig.  4  The  Q,  -  and  Q2  -modes  at  times  f3  =  7.5K  and  t4  =  20 K.  The  sharp  structures  in  the  Q2  -mode  behave  like  quasi- 
solitons  in  their  interaction  :  the  rapidly  moving  sharp  small  amplitude  structures  (with  varying  group  velocities  in  both  directions) 
pass  through  each  other  as  well  as  through  the  slowly  moving  large  scale  structures. 


Proc.  ofSPIEVol.  5105  279 


these  fine-scale  structures  are  transparent  not  only  to  other  fine-scale  structures,  but  also  transparent  to  the  slowly 
moving  sharp  large-scale  structures.  In  essence,  one  can  refer  to  this  as  soliton  turbulence.  The  coupling  to  the 
<2,  -mode  leads  to  the  Q,  quasi-solitons  being  broaden  further  with  the  emission  of  broadband  ’noise’.  The<22-large 
scale  structure  maxima  rapidly  pulsates  in  a  way  very  reminiscent  of  the  periodic  pulsating  bound  2-soliton  solution  for 
the  simple  NLS  [19]. 

In  Fig.  5  we  plot  the  <2,  -  and  Q2  -  modes  at  t  =  50  K  and  t  =  110  KL  to  contrast  their  different  time  evolution.  It  is 
apparent  that  the  group  velocity  of  the  large  scale  structures  in  both  modes  are  unchanged  by  the  coupling.  Both  modes 
exhibit  soliton  turbulence  while  the  <2i  -mode’s  major  2  structures  still  retain  their  identity,  albeit  broadened.  The  initial 
Q2  -mode’s  major  2  structures  (see  Fig.  2)  have  been  destroyed,  as  can  be  seen  in  Fig.  5  at  t  =  1 10  K. 
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Fig.  5  The  Qx  —  and  Q2  —modes  at  t  =  50 K  and  t  =  1 10 K.  The  major  peaks  in  the  Q2  —mode  pulsate  not  unlike  the  2-soliton 
bound  solution  for  simple  NLS,  except  that  for  the  Q2  —mode  the  major  peaks  have  a  non-zero  drift  velocity.  The  lower  amplitude 
peaks  that  are  superimposed  on  the  broadband  noise  act  like  soliton  turbulence. 


4.  CONCLUSION 

A  quantum  lattice  gas  representation  is  developed  for  the  solution  of  coupled-NLS  equations,  equations  that  can  model 
the  propagation  of  pulses  down  a  birefringent  fiber.  The  quantum  algorithm  requires  4  qubits  per  spatial  node,  with  the 
unitary  collision  and  streaming  operators  coupling  qubit  pairs.  The  qubit  pairs  are  then  coupled  to  each  other  by  the 
introduction  of  appropriate  phase  changing  potentials  to  the  qubit  pair  wave  functions.  Quantum  phase  coherence  needs 
to  be  enforced  for  only  1 6  quantum  gate  operations  since  measurements  need  then  be  taken  to  determine  the  absolute 
value  of  the  wave  functions  so  as  to  be  able  to  initialize  the  quantum  computer  for  the  next  time  evolution  stage.  It  is 
interesting  to  note  that  a  liquid  NMR  quantum  computer  [13]  has  recently  been  used  to  accurately  solve  the  diffusion 
equation  using  a  quantum  lattice  gas  algorithm.  Since  the  quantum  mechanical  sites  communicate  with  each  other 
classically,  this  architecture  has  been  dubbed  a  Type-II  quantum  computer  [6].  The  Type-II  quantum  computer  networks 
many  small  quantum  processors.  Given  this  preliminary  experimental  success,  it  appears  at  least  possible  that  the 
quantum  algorithm  presented  here  may  be  run  on  an  actual  Type-II  quantum  computer  in  the  future.  However,  the 
quantum  algorithm  for  vector  solitons  would  be  very  difficult  to  implement  experimentally.  Even  though  a  one¬ 
dimensional  model,  it  is  sufficiently  complex  because  amplitudes  must  be  extracted  from  the  quantum  measurement 
process  and  not  simply  probabilities  which  are  more  naturally  extracted  from  the  projective  von  Neumann  quantum 
measurement  process  over  a  large  ensemble  of  molecules  as  occurs  in  NMR  spectroscopy.  The  fact  that  amplitudes, 
instead  of  binary  bits  or  probabilities,  must  be  extracted  from  the  quantum  measurement  is  what  makes  the  experimental 
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implementation  much  less  efficient  than  otherwise.  Nevertheless,  it  is  hoped  that  this  one-dimensional  quantum 
algorithm  serves  as  a  stepping  stone  towards  more  complicated  quantum  algorithms  with  sufficiently  many  qubits  per 
site  so  that  the  experimental  trade-off  between  the  difficulty  of  achieving  quantum  entanglement  and  the  severe 
limitation  of  measurable  classical  information  that  can  be  extracted  from  the  quantum  computer  begins  to  work  in  our 
favor.  Finally,  it  is  interesting  to  note  that  the  quantum  lattice  gas  representation  is  very  efficiently  parallelized  on  a 
classical  computer. 
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